Quantum properties of counter-propagating two-photon states generated in a planar 

waveguide. 
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A nonlinear planar waveguide pumped by a beam orthogonal to its surface may serve as a versatile 
source of photon pairs. Changing pump-pulse duration, pump-beam transverse width, and angular 
decomposition of pump-beam frequencies characteristics of a photon pair including spectral widths 
of signal and idler fields, their time durations as well as degree of entanglement of two fields can be 
changed significantly. Using the measured spectral widths of the down-converted fields and width 
of a coincidence-count dip in a Hong-Ou-Mandel interferometer entropy of entanglement can be 
determined. 

PACS numbers: 42.50.Dv Quantum optics, 42.65.Wi Nonlinear waveguides, 42.65.Lm Parametric down- 
conversion and production of entangled photons 



I. INTRODUCTION 

The process of spontaneous parametric down- 
conversion as a source of entangled photon pairs has been 
used in numerous experiments during the last twenty 
years. The physicists went through a long way from the 
first pioneering experiments showing basic properties of 
entangled photon pairs [spectral (temporal) and polariza- 
tion correlations of photons comprising a pair, see, e.g. 
[U, 0]] to the recent sophisticated experimental setups 
that demonstrate quantum teleportation quantum 
cloning test Bell and other nonclassical inequalities 
[H'j'g'I, or generate Greenberger-Horne-Zeilinger states 7]. 
Usefulness of the fragile entangled photon pairs has also 
been verified in applications like quantum cryptography 
Q and absolute metrology ^ to name a few. 

Also sources of photon pairs have been improved sig- 
nificantly. Various geometric configurations of usual non- 
linear crystals (see, e.g., in [13, El) are gradually being 
replaced by new more efhcient sources based on quasi- 
phase-matching [H, [l^ [13, [Hi- Nonlinear photonic- 
band-gap fibers seem to be extraordinarily interesting as 
sources of photon pairs (emerging in the process of four- 
wave mixing) due to a high effective nonlinearity [l6|, [l7j . 
Nonlinear layered structures as sources of photon pairs 
are under investigation at present [H, [ll| . 

A great deal of attention has been devoted to the gen- 
eration of two-photon states with specific spectral prop- 
erties in bulk materials. Entangled two-photon states 
with coincident frequencies have been obtained using ex- 
tended phase-matching conditions, i.e. assuming group- 
velocity phase matching on the top of the usual phase 
matching [2lj, Such states have perfect visibil- 
ity in the usual interferometric setups and moreover al- 
low a very precise clock synchronization ^2^. Spectrally 
uncorrelated two-photon states have also been studied 
extensively because they seem to be extraordinarily use- 
ful for linear quantum computation that needs indistin- 
guishable photons with a perfect time synchronization 



[24 |25| . Although such states can be generated from 
usual bulk crystals for suitable crystal length and pump- 
beam waist in non-coUinear configurations [2^ , more 
flexible approaches have been suggested in '241 exploiting 
phase-matching in the transverse plane. Both coincident- 
frequency entangled and unentangled two-photon states 
can be obtained in a nonlinear crystal with achromatic 
phase matching, i.e. when pump-beam frequencies are 
decomposed such that every frequency propagates along 
a slightly different angle [1^, [1^ ■ Also the so-called non- 
linear crystal superlattices, i.e. structures composed of 
several identical pieces of nonlinear material and spacers, 
have provided additional degrees of freedom for tailoring 
properties of the generated two-photon states [s^, [Sll . 
Results appropriate for photonic-band gap structures |32l | 
are reached for a higher number of nonlinear pieces. 

Also planar nonlinear structures are promising as 
sources of entangled photon pairs because, using spe- 
cific geometric configurations, they provide large possi- 
bilities for tailoring properties of the generated photon 
pairs. One of the most p ersp ective configurations (sug- 
gested and elaborated in f33|, [13, [H, [s^ ) is based upon 
pumping a planar waveguide by a beam perpendicular 
to its surface. Signal and idler photons then emerge as 
counter-propagating guided waves. Disadvantage of this 
configuration is that the pump beam propagates through 
a very thin nonlinear medium which thickness is given 
by the depth of the waveguide. In order to suppress de- 
structive interference in the three- wave nonlinear process 
this thickness has to be of the order of pump-field wave- 
length. Thus, very low generation rates have to be ex- 
pected. To cope with a low efficiency of the nonlinear pro- 
cess more sophisticated structures have been suggested 
^37„ .38,] . They use Bragg mirrors both above and below 
the waveguide. Their properties are chosen such that 
the pump-beam electric-field amplitude is maximally en- 
hanced inside the nonlinear waveguide. This may re- 
sult in enhancement of the efficiency by several orders 
of magnitude. The first experimental demonstration of 
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this source has been already reported in (39|. The use 
of pump pulses with frequencies propagating along dif- 
ferent angles brings even more flexibility and so photon 
pairs with an arbitrary shape of a two-photon spectral 
amplitude can be generated [4C|. 

As shown in this paper using a simple model of pla- 
nar waveguide with parabolic index of refraction fssi , l4l| , 
properties of photon pairs (spectral and temporal widths 
of the down-converted fields and entanglement) gener- 
ated from this type of geometry can be modified in broad 
ranges simply by changing parameters of the pump beam 
(pump-pulse duration, pump-beam transverse width, an- 
gular decomposition of pump-beam frequencies). Spec- 
tral widths of the down-converted fields can spread from 
circa 1 nm up to several tens of nm. Spectrally un- 
correlated (separable) states as well as strongly entan- 
gled states can be observed. A method for the deter- 
mination of entropy of entanglement from the measured 
signal- and idler-field intensity spectra and width of the 
coincidence-count pattern in a Hong-Ou-Mandcl interfer- 
ometer is also suggested. 

The paper is organized as follows. Sec. II is devoted 
to the determination of a two-photon spectral amplitude 
of the generated photon pair. This amplitude is later 
used to derive spectral (Sec. Ill) and temporal (Sec. IV) 
properties of the down-converted fields, and characterize 
entanglement of the signal and idler fields (Sec. V). Ex- 
perimental determination of entropy of entanglement is 
discussed in Sec. VI. Conclusions are drawn in Sec. VII. 
Appendix A contains general formulas describing proper- 
ties of photon pairs. Appendix B is devoted to Schmidt 
decomposition of a two-photon spectral amplitude. 



II. TWO-PHOTON SPECTRAL AMPLITUDE 
OF A PHOTON PAIR GENERATED FROM A 
WAVEGUIDE 

We consider the generation of a photon-pair into 
guided modes of a planar waveguide made of LiNbOa 
that is pumped by a travelling-wave pump beam at the 
wavelength Ap — 1.064 x 10~^ m that propagates under 
the central angle 6^ with respect to the x axis orthogonal 
to the surface (see Fig. [T]). The pump beam is assumed 
not to be cross spectrally pure in general, i.e. differ- 
ent pump-beam frequencies can propagate under differ- 
ent propagation angles. The guided signal and idler fields 
then form counter-propagating beams. For simplicity, we 
pay attention to the waveguide with a parabolic profile of 
index of refraction n(x), n(x)'^ = 7i§(l — a^x'^) {a is pa- 
rameter of the waveguide) that supports only TE-guided 
modes with a gaussian profile. 

Energy of the nonlinear interaction that produces pho- 
ton pairs is described by Hamiltonian H: 



H{t) = eod / dV 




y-Z LiNbOs 



FIG. 1: Sketch of a nonlinear planar waveguide made of 
LiNbOs with parabolic profile of index of refraction along 
the X axis. The optical axis of LiNbOa as well as polariza- 
tion directions of signal, idler, and pump fields are parallel to 
the y axis. The signal (idler) field propagates along the +z 
( — z) axis with propagation constant jSs {Pi); the pump field 
with central wave-vector propagates under central angle 
with respect to the x axis. The waveguide is confined into 
the region {—Ly/2,Ly/2) in y direction. 



where E^"^ is positive-frequency part of the pump- 



beam electric-field amplitude and E^ '' {E\ ') stands 
for negative-frequency part of the signal- (idler-) beam 
electric-field amplitude operator. Symbol eo denotes per- 
mittivity of vacuum, d is effective second-order nonlinear 
coefficient, and h.c. means a hermitian conjugated term. 
Integration in Eq. ([1]) is over interaction volume V. 

In the considered waveguide, positive- frequency parts 
of the signal- and idler-beam electric-field amplitude op- 



erators Es and E^ ' can be decomposed as [Ep 
(4+V]: 

i?^+^(r,i) = dujaeair,uja)aa{uja), 



(-) _ 



eair,UJa) = C'a(Wa)reCt(_L^/2,L„/2)(2/) 



X exp 



\Ca{0Ja)\' 
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exp(±i/3Qz) exp(-iwaOi 

a = s,i; (2) 



symbol aa(wa) denotes annihilation operator of a mode 
with frequency LOa in field a. The sign -I- (— ) in the sec- 
ond relation in Eq. ^ is for the signal (idler) field that 
propagates along the (— z) axis. In Eq. ([2]), no(ti'a) is 
index of refraction of field a with frequency Wo, c means 
speed of light in vacuum, h, reduced Planck constant, Ly 
width of the waveguide along the y axis where a rectan- 
gular profile is assumed, and 7a (wa) = \/ n[){uja)LOaCi/ c- 
Function rect(o^b) {x) equals 1 for a < x < b and is zero 
otherwise. Normalization constant Ca in Eq. ^ has been 
determined from the condition that a photon emitted into 
a guided mode has energy huja' 

2eonl{uJa) dV\ea{r,uJa)\''^ ^ huja, a^s,i. (3) 



(1) We assume that the interaction volume V has 'length' 
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c/no(ajJ^) along the z axis when determining constant Cq 
in Eq. ©. 

Propagation constant Pa of field a along the z axis is 
given as follows (for details, see (41i]): 



ac 



C V nQ{uJa)uJa 

It can be approximately expressed as: 



(4) 



(5) 



where Va denotes group velocity of field a. 

The travelling-wave pump beam with central frequency 
and propagating along central angle 6^ is assumed to 
have a gaussian profile along the y and z axes character- 
ized by widths Yp and Zp, respectively, and is also gaus- 
sian in the time domain with pump-pulse duration Tp 
and chirp parameter a^. Different monochromatic com- 
ponents of the pump beam can propagate along different 
angles 9p{ujp), e.g. as a consequence of pump-beam reflec- 
tion on an optical grating or after propagation through a 
prism. The pump-beam positive-frequency electric-field 
amplitude E'p^^ can be written in the form: 



dujpE^'^\r,ujp) exp{-iujpt), 



Y2 

p 



El+\r,ujp) = Cp{ujp)exp (^-^^ exp 

X exp[ikp sm{dp{LUp))z] exp[— ifcp cos{9p{u!p))x] 
Tl{up-ulY 



X exp 



4(l-fmp) 



2Tnreonl{ujp)YpZp{l + aj) t;^ cos(6'p(wp))/ 

(6) 

Vp is pump- field group velocity [kp{LLip) = + (ujp — 
^p)/vp, fcp = fcp(^p), = dkp/dujp\^^^^o], Pp pump- 
field power, and / denotes repetition rate of the pulsed 
pump field. 

First-order perturbation solution of the Schrodinger 
equation using Hamiltonian H{t) given in Eq. ([T]) and 
assuming the incident signal and idler fields in vacuum 
states provides the following expression for an outgoing 
two-photon state j?/;'^^): 

IV'^^^) = / dujs / duJ,<^>^P{LJs,u^)al{us)d\{uJi)\va.c). (7) 



Two-photon spectral amplitude ^^^(tJs, w,) giving the 
probability amplitude of having a signal photon at fre- 
quency ljJs and an idler photon at frequency uji generated 



from one pump pulse is derived in the form: 
.^/2^T2^T'^eQd 



UJ. 



Cp{uJs+OJi)C*{uJs)C*{uj{) 



_ erf exp 



4(1 -I- iap) 



X exp 



X exp 



Zl[kp siTL{9p{ujs + uji)) - Ps ^ Pi 



k'^cos'^{ep{ujs +tOi)) 



(8) 



symbol erf stands for error function [erf(a::) = 
2/V^/o'exp(-y2)dy]. 

Using first-order Taylor expansions for propagation 
constants kp, /3s, and f3i as well as for propagation an- 
gle Op together with second-order Taylor expansion for 
the expression 1/(72 -|- ^f) [ 1/(72 + jf) w go -f guAuj^ + 

a = s, i] we arrive at a two-photon spectral amplitude $ 
giving contribution from / pump pulses and having the 
following gaussian form: 



ZjiTr^ 



^{u)s,uji) = C$ J— ^-^ exp[-0(a;s,t^»)]. 

(/)(ws, Wi) = !islA^1 + IiilA^l + /2si 

+ /isAc^, + /i,Ac^, + /o. (9) 

Coefficients / occurring in Eq. ([9]) are expressed as fol- 
lows: 
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z'- 
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Qa 



4(1 -hwp) ' 4 
_ fcOcos(^O)/ 



klcOs{el)g2a 

klsin{el)be^ 



gia 



cos(0P) 4sin(^g) ^^ 
k^v?: Vr, ^ 



fcgcos(2g°) 
cos(0O) 



90 



f2si — 



VpsVpiZp 



2{1 + iap) 2 

fc>s(^/ 
Ssi = [ Kcos{ei)g2. 



+ Qsi, 



COS(0O) 



klsin{6l)be^ 



(gis + gii) 



cos{e°) 4sin(6lO)- 
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fc°cos(2g°) 
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/la 



fc^cos(6'P ^ ^ gia 



COS 



50 



/o 



[fcOcos(^O)]^ 



■50, 



(10) 



and 



Vps = 



sin(0O) 



sin 



fc;^cos(0O)De,-- 



in(0O) 



Vn 



1 
1 

n,. 



(11) 



Coefficient describes angular decomposition of 

pump-beam frequencies; Dg^ — d0p{ujp) / dujp\^ ^^c,; 

Dg^ = I?ep(cjp)^/(27rc). Influence of frequency filters 
with a gaussian shape is described by their widths as 
and fTi (for the signal and idler fields, respectively) that 
occur in Eq. PU)) . 

Normalization constant C$ introduced in Eq. Q is 
determined along the expression: 



Pn 



(12) 



Phase matching for central frequencies has been as- 
sumed when deriving the expression for two-photon spec- 
tral amplitude $ written in Eq. ([9]), i.e. 



fc° sin(0°) A" =0. 



(13) 



Equation (|T3l) represents condition for possible values of 
central frequencies uj^, a;° and central angle 6^ of 
pump-beam propagation. This condition even with the 
inclusion of quasi-phase matching has been extensively 
studied in ^55*1. 

The role of pump-pulse duration Tp and pump-beam 
transverse width Zp on the shape of two-photon spectral 
amplitude $ can be understood when we transform the 
amplitude $ into new variables fl and ui; ^ {uJs+ij->i)/2, 
u; = {lOs - uJi)/2: 



^in,u;) = 2C, I ^P""" 



X exp'^ 



l + «p 



1 -I- iap 



Zl{Vps + Vpi)Vs^ 



ZliVps + Vp^f 
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(14) 



(15) 



(16) 



Coefficients Qs, Gi, and Gsi occurring in Eq. (fTOj) have 
been neglected when the expression in Eq. has been 
derived because they are small in comparison with those 
written explicitly in Eq. (jl4|) under the studied condi- 
tions. We can see from Eq. that pump-pulse dura- 
tion Tp influences only the coefficient of quadratic form in 
sum frequency Q, whereas pump-beam transverse width 
Zp and widths of spectral filters as and at modify all of 
them. The shape of a two-photon amplitude can also 
be controlled using parameter Dg of angular decompo- 
sition of pump-beam frequencies that occurs in expres- 
sions for the coefficients multiplying Afi^ and AOAw in 
Eq. (fH|) . Nonzero values of parameter Dg^ lead to rota- 
tion of the shape of the two-photon amplitude $ in the 
plane spanned by frequencies u/g and uji 2&., .29] . 

Considering frequency degenerate case (w° = uj°, i.e. 
and Op =0), cross spectrally pure pump beam 
0), and omitting frequency filters {as, at —^ oo) 



Vs = v., 

{Dg — U), and omitting trequency niters (as , a 
we obtain the two-photon spectral amplitude $ in a sim 
pie form: 



^{n,uj) 




X exp 



A17" 



ZW^ 



;(i7) 



i.e. pump-pulse duration Tp determines properties de- 
pending on sum frequency Q whereas pump-beam trans- 
verse width Zp is responsible for properties related to 
difference frequency lu. This behavior is similar to that 
occurring in coincident-frequency entangled two-photon 
states generated from bulk materials and described, e.g., 
in [20j (length L of a crystal plays the role of Zp). 



III. SPECTRAL PROPERTIES OF PHOTON 
PAIRS, PHOTON-PAIR GENERATION RATE 

Number N of photon pairs generated in 1 s is deter- 
mined along the formula: 



N = duOs I duj^\^{uJs,^^i)Y 



(18) 



The general expression for the number N of photon pairs 
assuming the spectral two-photon amplitude $ in the 
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form written in Eq. ([9]) can be found in Appendix A 
[Eq. (|Aip ]. Neglecting coefficients Gs, Gi, and Gsi in 
Eq. (fTO|l . a simplified expression can be derived: 



(19) 



where 



4(1 + a2) ^ 







) {<^'s 




p *^ pz 
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(20) 



If frequency filters are omitted, the expression in Eq. JT 
for photon-pair generation rate N gets a simple form: 



1 + alVs, 



i.e. the generation rate N does not depend both 
on pump-pulse duration Tp and pump-beam transverse 
width Zp. This is a consequence of geometry of the 
considered three-mode interaction. Assuming values of 
the waveguide parameters as written in Fig. [2] and inci- 
dent pump-field power Pp = 1 W, the generation rate 
N equals 3 x 10^ s^^, i.e. if the pulsed pumping has 
repetition rate f — 8 x 10^ s^^, a photon pair is gener- 
ated from one pump pulse with probability 3.8 x 10"'*. 
Taking into account the depth of the waveguide of the 
order of pump-field wavelength, this probability is high. 
There are two reasons. The nonlinear process exploits 
the largest element of the nonlinear tensor d of LiNbOa 
(this gives 2 orders of magnitude in comparison with com- 
monly used orientations). Also the down-converted fields 
are confined in their transverse profiles into very narrow 
regions, so the electric-field amplitude per one photon is 
high which improves efficiency of the nonlinear process. 

Signal-field intensity spectrum Ss determined accord- 
ing to the formula 



(22) 



takes a gaussian form: 



Ssic^s) = Ss exp 



(23) 



Amplitude Sg , width a^j^ , and shift Suj'^ of the center 
for the signal-field intensity spectrum assuming negligible 
coefficients Gs, Gi, and Gsi are given as: 



Ss = |C*pexp(-2/o) 



1 + al 



2(l + a^)+af 



Z^V^. 

P * pi 



-, -1/2 



(24) 



2(1 + a2) 



2 ' 



(25) 
(26) 



DfT is given in Eq. (|20p . The general expressions for 
parameters Ss, a^i, , and Saj^ can be found in Appendix A 
[Eqs. (|A4| — 1X6]) ]. Characteristics of a gaussian idler-field 
intensity spectrum Si{uJi) can be derived from symmetry. 
When frequency filters are not included, the expression 
in Eq. for width (7^^^ can be further simplified: 




(27) 



This means that the signal-field spectral width a^^ (and 
similarly the idler- field spectral width aoj.) decreases 
with increasing pump-pulse duration Tp and pump-beam 
(21) transverse width Zp. Assuming cw pumping, spectral 
widths and a™ are inversely proportional to pump- 



beam transverse width Zr, 



V2 

VsiZp 



(28) 



On the other hand the following expressions hold for the 
pump-beam transverse width Zp sufficiently large, 



V^\VpA^ 



1 -h ftp 



I Zr, — *00 



Vs,. 



V2\Vps 



l + a2 



Vs, 



(29) 



i.e. widths a^^^ and a^^^ are inversely proportional to 
pump-pulse duration Tp. Dependence of the width of 
the signal-field intensity spectrum \ax^ = 2-kc/ {uu^sYa^^] 
on pump-pulse duration Tp and pump-beam transverse 
width Zp is shown in Fig. [2]for the considered waveguide. 
Narrow spectra broad circa 1 nm are generated if pump- 
pulse duration Tp is sufficiently long and also pump-beam 
transverse width Zp sufficiently wide. On the other hand, 
intensity spectra wide several tens of nm can be observed 
for ultrashort pump-pulses and strongly focused pump 
beams. The reason for this behavior is that shortening 
of a pump pulse and focusing of a pump beam lead to 
weakening of frequency and phase matching conditions. 

Ratio (Juisl^i^i of spectral widths of the signal- and 
idler-field intensities is important in some applications 
[2^ . For example, photon pairs used in heralded single- 
photon sources should preferably be composed of one 
photon with a narrow spectrum (convenient in propa- 
gation through an optical fiber) and one photon with 
a wide spectrum (leading to high detection efficiencies 
when post-selecting). Provided that the pump beam is 
cross spectrally pure, the ratio o^^Jo^^^ is given by mate- 
rial constants (group velocities of three fields) for given 
values of pump-pulse duration Tp and pump-beam trans- 
verse width Zp. The ratio a^^/a^. equals 1 if the signal 
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FIG. 2: Contour plot of width of the signal-field intensity 
spectrum as a function of pump-pulse duration Tp and pump- 
beam transverse width Zp is shown; logarithmic scale is used 



on the z axis; X° = 1.064 x 10""" m, A' 
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0, K 



1 X lO"-" m, Ly 



1 X 10" 



Dflp = deg m-\ a = 4 x lO" 
41.05 x 10"" mV"^ 



cx), d — 



and idler fields are symmetric (w" ~ ujf). On the other 
hand, parameter Dg characterizing angular decomposi- 
tion of pump-beam frequencies can substantially change 
this ratio. Either narrowing or broadening of an inten- 
sity spectrum of a down-converted field can be reached 
by the change of values of parameter Dg^ , as formulas in 
Eqs. ((27)) and pT|) indicate. Values around 10 for the 
ratio <yujs/<^LJi can be obtained for the waveguide with 
values of parameters defined in Fig. [21 as documented in 
Fig. [21 We note that high values of this ratio occur when 
the two-photon spectral amplitude <I> is nearly spectrally 
uncorrelated (compare Fig. [51 later). 



IV. TEMPORAL PROPERTIES OF PHOTON 
PAIRS 

Temporal properties of photon pairs can be conve- 
niently described using two-photon amplitude $ in the 
time domain that is given as a Fourier transform of that 
in the frequency domain: 



1 

2^ 



x exp(— zcJsTs) exp(— jWiTi). (30) 



Because the two-photon spectral amplitude <i> as given in 
Eq. ([9]) is gaussian, the two-photon amplitude $ defined 
in Eq. ([50)1 takes also a gaussian form, that can be found 
in Appendix A [Eq. ()A10[) ] . 

Photon flux Ns in the signal field is then determined 
along the formula [l^ 



(31) 



FIG. 3: Contour plot of ratio a\J(T\^ of spectral widths of 
signal- and idler-field intensities as a function of parameter 
Dg"' giving angular decomposition of pump-beam frequencies 
outside the waveguide and pump-beam transverse width 
Zp is shown. Ratio (J\^/(J\^ is plotter for Dg"' < 0. The 
bold curve indicating spectrally uncorrelated states is given 
using the formula in Eq. ([55)l later. We have D^"* = 

sin(^p)/cos(^°"*) dnoiujp) /dujp\^^^^a and 

arcsin[no(ajp) sin(Sp)]; Tp = 1 x 10" 
other parameters are given in Fig. [51 



gout 



s, values of the 



and attains the following gaussian form; 



Ns{ts) = Us exp 



(32) 



Considering a pump pulse without chirp (op — 0) and 
omitting contributions in Eq. (|10|) given by coefficients 
Qs^ Gi, and Qsi, we arrive at the following simplified ex- 
pressions: 



5r^ 




1 -1/2 



(33) 

(34) 
(35) 



coefficient 2?/ is given in Eq. ()Alip in Appendix A. The 
formula in Eq. ([34)) shows that the signal-field tempo- 
ral width ffr^ increases as the pump-pulse duration Tp 
and pump-beam transverse width Zp increase. Also the 
narrower the signal-field frequency filter, the greater the 
temporal width . 

Depending on pump-pulse duration Tp and pump-beam 
transverse width Zp width CTt-^ of the signal-field intensity 
profile spreads from several tens to several hundreds of fs 
(see Fig. [J) . Width ar, of the signal-field intensity tem- 
poral profile has to be greater than that characterizing 
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FIG. 4: Contour plot of temporal width ctts of signal-field 
photon flux as a function of pump-pulse duration Tp and 
pump-beam transverse width Zp is shown; values of the used 
parameters are given in Fig. O 



the pump field (rp/ \/2) as the expression in Eq. con- 
firms. Also the greater the pump-beam transverse width 
Zp the greater the width as a consequence of propa- 
gation of a signal photon along the waveguide. Angular 
decomposition of pump-beam frequencies {Dg ) can ei- 
ther broaden or shorten the pulsed photon fluxes Ng and 
Ni as follows from formulas in Eqs. ([M]) and pT|) . 

Temporal and spectral widths of the down-converted 
fields are not independent provided that the pump pulse 
is not chirped (op = 0): 



(36) 



Considering the symmetric case (w^ — lo^), cross spec- 



trally pure pumping [Dg 
and no frequency filters [us 



0) without chirp 



0), 



fJi — !■ oo), the following rela- 



tion for the product of temporal and spectral widths of 
the down-converted fields can be derived: 



1 

Ct = — 

• ^" 2 



-'a ' p 



> 1, 



(37) 



Product (Tctj^cTTa written in Eq. p7p goes to infinity for 
cw pumping. 

Temporal properties as described by photon fluxes Ns 
and Ni are not accessible directly experimentally owing 
to very short time scales. Thus interferometric setups are 
needed to obtain these characteristics. Temporal overlap 
of the signal- and idler-photon wave- functions can be de- 
tected in a Hong-Ou-Mandel interferometer that provides 
a normalized coincidence-count rate as a function of 
mutual time delay rj between the signal and idler photons 
El: 



R^Ti) = 1 - pin), 



(38) 



where 



1 /■°° 

PiTl) = JT dtA 



dtf 



Re{^{tA.tB-Ti)<^*{tB,tA-n)} (39) 

and the number N of photons generated in 1 s is given by 
Eq. p^ . Symbol Re means the real part of an argument. 
The shape of two-photon amplitude $ thus determines 
the pattern of coincidence-count rate i?„(T;) as discussed, 
e.g., in [H. 

Using the two-photon spectral amplitude <f> written in 
Eq. ([9]) the normalized coincidence-count rate i?n can be 
derived in the form: 



i?„(rO 1 - Aexp(-BTf) cos[(^° - L^°)r,] 



(40) 



Coefficients A and B arc given as follows neglecting co- 
efficients Qg, Qi, and Qsi [the general expressions can be 
found in Eqs. (|A2ip and (|A22|) in Appendix A]: 



Zl{l + al) (K 



ps 



*plJ 



4r2 



B = 



^p ^ SI 



y2 

SI 

-i -1 



-1/2 



(,41) 
(42) 



Visibility of a coincidence-count pattern described 
by the formula in Eq. ((40|) and defined along the expres- 
sion (i?n,min givcs a minimum value of Rn) 



V 



Rnin ^ oo) - i?n 



is obtained in the form: 
V 



A 



2- A 



(43) 



(44) 



For cw pumping (rp — > oo) and without frequency filters, 
coefficient A — > 1 and also visibility V ^ 1. Visibility V 
equal to 1 is also reached for the symmetric case (w^ — w,^ 
leads to A = 1). 

On the other hand, coefficient B given in Eq. (|^^ de- 
termines width Ati of the coincidence-count dip defined 
by the condition i?n(Ar;/2) = l — A/2. This relation can 
be rewritten into a transcendent equation, 



exp 



-An 



0, 



(45) 



that has solution in the range An G (0,27r/[a;° — a;°]). 
We note that oscillations at frequency uj^ — w° occur in 
the coincidence-count pattern of rate Rn- Coefficient B 
depends only on pump-beam transverse width Zp and 
widths as and ai of frequency filters. The narrower the 
frequency filters and the wider the pump-beam transverse 
width Zp, the smaller the value of coefficient B and also 
the larger the width At; of normalized coincidence-count 
rate i?„. We can see that width An of the coincidence- 
count dip does not depend on pump-pulse duration Tp 
and pump-beam parameter Dg . The reason is that a 
Hong-Ou-Mandel interferometer monitors only the differ- 
ence of the signal- and idler-field frequencies 2uj = ujs — oJi 
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200 n 




I 

5 25 45 65 85 105 



Zp (10"'m) 

FIG. 5: Width Ar; of the dip in normaUzed coincidence-count 
rate Rn in a Hong-Ou-Mandel interferometer as it depends on 
pump-beam transverse width Zp; values of the used parame- 
ters are written in Fig. (2] 



symbol ln2 means logarithm of base 2. Using Eq. (|47|) we 
arrive at: 

^. = -log,(l-.)-i^. (49) 

As for possible values of entropy Se, there are two 
boundary cases. If d ^ 1 (also P and |e2| = e2c) 
then all eigenvalues A„ are equal, i.e. we have a max- 
imally entangled state. On the other hand, if 
(also P oo and e2c 0) there is only one nonzero 
eigenvalue Aq; i.e. the two-photon spectral amplitude 
$(tJs7'^i) factorizes and describes a separable state use- 
ful, e.g., in linear quantum computation [23|. 

A. Entangled states 



that depends only on pump-beam transverse width Zp in 
the considered configuration [compare Eq. (fH)) ]. 

Width Ati of the coincidence-count dip can be con- 
trolled changing pump-beam transverse width Zp in a 
broad range of its values, as documented in Fig. O This 
property may be useful in metrology applications (e.g., in 
measurement of modal dispersion of a waveguiding struc- 
ture). 



V. ENTANGLED AND SEPARABLE STATES 

Entanglement between the signal- and idler-field fre- 
quencies can be conveniently quantified using entropy of 
entanglement [H, 13]. To determine entropy Se of en- 
tanglement we have to decompose a two-photon spectral 
amplitude <f> into Schmidt decomposition p^ : 

^(WsjWi) = ^Xn4>s,n{<^s)<Pi,niuJi), (46) 
n 

where A„ are coefficients of the decomposition and func- 
tions 0s,Ti(ws) a-nd (pi^ni^i) form the Schmidt basis. This 
decomposition assuming a gaussian two-photon spectral 
amplitude $ is found in Appendix B and gives coefficients 
A„ in the form of geometric progression: 

A„ = n = 0,l,...,c5o. (47) 

Value of parameter § is derived from values of param- 
eters characterizing amplitude $ using Eqs. (jB8[) and 
(|B11|) in Appendix B. We note that Yl,n -^n = 1 as a 
consequence of assumed normalization of the amplitude 
$; / dujs^ dhJi\^{ujs,u)i)\'^ = 1. 

Entanglement between the signal- and idler-field fre- 
quencies can be determined using entropy Se derived 
from eigenvalues A„ of the Schmidt decomposition, 

Se^~Y.^l\og^{\l); (48) 

n 



Entangled states with high values of entropy Se of en- 
tanglement are generated if either pump-pulse duration 
Tp is short or pump-beam transverse width Zp is narrow, 
as the analysis contained in Appendix B shows. The 
shorter the pump-pulse duration Tp the greater the value 
of entropy Se- The narrower the pump-beam transverse 
width Zp the greater the value of entropy Se- Also the 
wider the widths tJs and ai of frequency filters the greater 
the values of entropy Se- 

On the other hand, widths a^}^ and au}^ of the signal- 
and idler-field intensity spectra increase with decreas- 
ing pump-pulse duration Tp and pump-beam transverse 
width Zp- The used frequency filters make these spectra 
narrower. Comparing qualitatively this behavior with 
that of entropy Se of entanglement described above, we 
can conclude that the wider the spectra of the signal and 
idler fields, the better the entanglement of the signal and 
idler fields. 

Instead of characterizing entanglement by entropy Se 
we can judge it by a minimum number nmin of eigen- 
functions (modes) from the Schmidt basis that restore a 
two-photon spectral amplitude $ with probability Pmin- 
The number nmin satisfies the following inequalities: 

\l < Pmin A ^ A^ > p^i„. (50) 

11=0 n=0 

The value of probability Pmin should be chosen with re- 
spect to the precision of measurement. 

For the considered waveguide, entropy Se of entan- 
glement as well as minimum number rimin of modes are 
shown in Fig. [S] as they depend on pump-pulse duration 
Tp and pump-beam transverse width Zp. We can see in 
Fig. [6] that entanglement of the signal and idler fields is 
rather weak in a broad area around the bold curve char- 
acterizing spectrally uncorrelated (separable) states and 
so we can approximate the generated state by a separable 
two-photon spectral amplitude Entangled states for 
which several modes in the Schmidt decomposition are 
necessary occur on the borders of the contour plot, i.e. 
where the quantity TpV^/Zp considerably differs from one 
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a) 




b) 




Z (10 m) 



FIG. 6: Contour plots of entropy Se of entanglement (a) and 
minimum number nmin of modes in Schmidt decomposition 
(pmin ~ 0.95) (b) as they depend on pump-pulse duration 
Tp and pump-beam transverse width Zp; bold curves in the 
middle of the plots are given by the condition in Eq. (|53|l for 
spectrally uncorrelated states; values of the used parameters 
are given in Fig. [2l 



[see Eq. (ISBl) later]. Wide intensity spectra occur in this 
region (compare Fig. [2]). 

The narrower the frequency filters the lower the val- 
ues of entropy and thus the weaker the entanglement 
between the down-converted fields, as demonstrated in 
Fig. [71 Assuming Pmin — 0.95 and values of the waveg- 
uide parameters appropriate for Fig. [3 filters narrower 
than 12 nm transform the down-converted fields into a 
separable state whereas two independent modes are suf- 
ficient for wider filters. 



B. Spectrally uncorrelated (separable) states 

The condition for a separable state derived from the 
formula in Eq. e2c — 0, is fulfilled provided that 

/2si = as follows from the definition of coefficient e2c in 
Eq. (jB4p in Appendix B. Separability of the two-photon 
spectral amplitude $ written in Eq. ^ is clearly visible 
in this case and shows that separable states are generated 



CO 



1.0 

0.8 
0.6 
0.4 



0.2 z 



0.0 




a (10"V) 

FIG. 7: Entropy Se of entanglement as a function of width a 
{as = ai = a) of frequency filters in nm; Tp = 1 x 10^^^ s, 

: 5 X 10"'' 
eters are written in Fig. [2l 



only provided that the axes in which the quadratic form 
given in Eq. is diagonal coincide with the ojs and tUi 
axes. Eigenvalues fii^2 of the quadratic form written in 
the second line of Eq. ^ are given as follows: 



Ml, 2 



.f2s + ^2^ ± V(/2s - f2^Y 



2si 



(51) 



whereas angle giving declination of the axes of the di- 
agonal quadratic form from the and LOi axes is written 
as: 



tan(V...) = /^--/^'^ + (52) 



2si 



The limit f2si in Eq. ((5^ leads to ipsi = 0, i.e. the 
axes of diagonal quadratic form coincide with the uj^ and 
LOi axes. 

Substituting expressions in Eqs. Uni) and ^ into the 
separability condition J2si = we arrive at: 



klcos{Bl)be 



1 



2gs^ = 0. (53) 



Assuming fixed values for pump-pulse duration Tp and 
pump-beam transverse width Zp and no chirp (op = 0) 
the condition in Eq. (|53p represents a quadratic equation 
for parameter Dg^ of angular decomposition of pump- 
beam frequencies and its solution takes the form: 



1 



1,2 



2A;Ocos(0O) 



2sin(^0) ^ ^ 




(54) 



Solution for Dg written in Eq. (I54p exists only when 
argument of the square root in Eq. (|54p is non-negative. 
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(deg) 




FIG. 8: Angle ij^si giving orientation of axes of the diago- 
nal quadratic form of two-photon spectral amplitude $ with 
respect to the ujs and oJi axes as it depends on pump-beam 
transverse width Zp and parameter -D^"* giving angular de- 
composition of pump-beam frequencies outside the waveguide 
(for its determination, see Fig. [Sjl . The bold curve indicates 
spectrally uncorrelated states and is given by the formula in 



Eq. (|55|) : Tp = 1 x 10 s, values of the other parameters are 
the same as in Fig. (2] 



VI. EXPERIMENTAL DETERMINATION OF 
ENTROPY OF ENTANGLEMENT 

If the pump field is not chirped (op = 0), a gener- 
ated photon pair is described by a two-photon spectral 
ampHtude <I> in a gaussian form written in Eq. ^ that, 
apart from central frequencies, is specified by three real 
parameters. These parameters can be experimentally de- 
termined measuring widths of the signal- {a^^^ ) and idler- 
field {<Jui ) intensity spectra and width At; of coincidence- 
count pattern in a Hong-Ou-Mandel interferometer, as 
shown below. 

Coefficient B describing width of the coincidence-count 
dip in the Hong-Ou-Mandel interferometer can be ob- 
tained by fitting the experimental curve Rtx{ti) using the 
prescription written in Eq. (j40p . The measured widths 
Gi^^ and a^^ of signal- and idler-field intensity spectra 
provide coefficient F [for definition, see Eq. (|A9p in Ap- 
pendix A]: 



F = 



(57) 



Thus, there is a lower limit for possible values of pump- 
beam transverse width Zp keeping pump-pulse duration 



Tp fixed. 



Considering the symmetric case (cj^ = ujf) and ne- 



glecting coefficient Qs 
simplifies: 



the formula in Eq. ([5^ for Dg 



± 



1,0 

ftp 



'p 



(55) 



and is valid provided that Zp > VsTp. If the pump beam 
is cross spectrally pure, i.e. Dg = 0, the following con- 
dition assures the generation of a separable state: 



Further, Eqs. (|A9l) and (|A22|1 in Appendix A can be re- 
cast into the form: 



fr 

J2i 

J2si 



fr 
J 2s 

F ' 

F + l 



F 



I 2s 



1 

2B' 



(58) 



Substitution of Eqs. ^ into Eqs. and (|X5D in Ap- 
pendix A leads to a quadratic equation for coefficient 



F 



(JLf 



Solution of Eq. ([55]) gives the value of coefficient /j^. 
the symmetric case (cj^ = lj^), F ^ \ and we have 



For 



1 

Vs 



(56) 



rr 

J2s — 



8B{a^2 - B) ■ 



(60) 



Changing values of parameters Tp, Zp, and Dg^ photon 
pairs with an arbitrary gaussian two-photon spectral am- 
plitude $ can be generated [J^] . An arbitrary orientation 
of axes of the diagonal quadratic form of two-photon am- 
plitude $ can be reached changing the value of parameter 
Dg^ of angular decomposition of pump-bcam frequencies, 
as demonstrated in Fig. [S] Parameters Tp and Zp then 
control spread of the two-photon spectral amplitude $ 
along the axes given by angles ipsi and ipsi + tt/2. Com- 
parison of graphs in Figs. [3] and [8] with special attention 
to frequency uncorrelated states leads to the conclusion 
that the greater the pump-beam transverse width Zp the 
greater the ratio cr^^/cr^; of intensity spectral widths of 
the down-converted fields. Such states have been found 
useful in quantum communication protocols. 



Values of coefficients f2i and f2si are then given by for- 
mulas in Eqs. (|58p . Coefficients /is and fn occurring in 
Eq. ^ give only shifts on the frequency axes lUs and uji 
and can be neglected. Also coefficient /o in Eq. ^ giving 
normalization can be omitted. 

Knowing values of coefficients /2s, f2i, and f2si, en- 
tropy of entanglement can finally be determined along 
the formulas given in Eqs. dH, (IbTT|1 . (IB8|1 . and ([49]). 



VII. CONCLUSIONS 

Properties of the down-converted fields can be con- 
trolled in wide ranges of values of characteristic param- 
eters using pump-pulse duration, pump-beam transverse 
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width, and angular decomposition of pump-beam fre- 
quencies in a waveguide with counter-propagating down- 
converted fields and transverse pumping. Widths of in- 
tensity spectra of the down-converted fields may vary 
from nanometers to tens of nanometers. Durations of 
the down-converted pulsed fields extend from tens of fs 
to several ps. Both entangled and separable (spectrally 
uncorrelated) photon pairs useful in linear quantum com- 
putation can be generated. Also attainable widths of a 
coincidence-count dip in a Hong-Ou-Mandel interferom- 
eter lie in a broad interval that is of interest in metrol- 
ogy applications. Using the measured spectral widths of 
the signal and idler fields and width of the coincidence- 
count dip, entropy of entanglement as well as parameters 
characterizing a two-photon spectral amplitude can be 
obtained. 

The considered nonlinear waveguide is promising as a 
versatile source of photon pairs in near future provided 
that efficiency of the nonlinear process is increased. 
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APPENDIX A: GENERAL FORMULAS FOR 
PHYSICAL QUANTITIES CHARACTERIZING 
THE EMITTED PHOTON PAIRS 

1. Spectral properties 

Substitution of the general form of two-photon spectral 
amplitude $ as given in Eq. ^ into Eq. results in 
the following expression for the number N of generated 
photon pairs in 1 s: 



|C$pexp(-2/o) 



nZpTp 



(l + a2)yp7 



where 



and 



£r £2 \ £r f2 £r £ £ 

J2sJli^ J2iJls J2siJ'^sJU 



(Al) 



(A2) 



exp 2 



(A3) 



Superscript r indicates the real part of a given complex 
coefficient. 

A gaussian signal- field intensity spectrum 6*^(^3) as 
written in Eq. (|23p is determined by its amplitude s^, 



width a^^ , and shift 5uo^^ of the center that are given in 
general by the following expressions: 

|C.pexp(-2/o)4^;^4=£/., (A4) 



V2(l -f al) 




(A5) 
(A6) 



where coefficients Vjt and are given in Eqs. (jA2p and 
(IA3p . Expressions for amplitude Si, width CTi^. , and shift 
Slu'^ belonging to the idler-field intensity spectrum Si{uji) 
can be obtained from those written in Eqs. (|A4| — [X6)) by 
an exchange of indices s and i. 

Relations among parameters of the signal- and idler- 
field intensity spectra can be established for a two-photon 
spectral amplitude <J> written in Eq. 



^? r- 
1 



Symbol F, 



F : 



F 



J 2s 
£r ' 
J2i 



(A7) 
(A8) 

(A9) 



gives the ratio of parameters that characterize a gaussian 
two-photon spectral amplitude 



2. Temporal properties 

Substituting the expression in Eq. Q for two-photon 
spectral amplitude ^{uJs,i-Ui) into the definition of ampli- 
tude $(rs,Ti) occurring in Eq. pO|) . we arrive at 



X exp 



X exp 




(1 + al)Vf 
ifisf + f2s{n - ifiif 



f2si{js ~ ifls){Ti ~ ifli) 



and 



Vf = Ahsh^ - fl 



(AlO) 



(All) 



Parameters of the gaussian form of signal-field photon 
flux Ng written in Eq. p2p [amplitude Us, width , and 
shift STg] are given for the general form of two-photon 
spectral amplitude $ in Eq. ([9]) by the following formulas, 
similarly as in the case of intensity spectra: 



= |C$pexp(-2to) 



V2{l + al) \Vf\Vhl 
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and 



Vt = 



2^21^15 — t2sitli 



exp 2 



2il-ls 



(A12) 
(A13) 

(A14) 

(A15) 
(A16) 



Coefficients t occurring in the above Eqs. (|A12I — [A16p 
are given as: 

t2s = Re{/2,/I?/}, 

t2r = Re{/2s/I?/}, 
t2si = -^e{f2si/Vf}, 

= Im{(2/2ji, - /2.Ji,)/I?/}, 
tu = Im{(2/2,/i,-/2.^/is)/2?/}, 

to = /o-Re{(/2,/L + /2s/l'.-/2«Jl«/H)/2?/}; 

(A17) 

Re (Im) denotes a real (imaginary) part of an expression. 
Similarly as in the case of spectral properties, relations 
among parameters of the signal and idler fields can be 
derived: 



to. 



Parameter T, 



T = — 

t2s 



(A18) 
(A19) 

(A20) 



gives the ratio of parameters that characterize a gaus- 
sian two-photon amplitude $ in the time domain. If the 
pulsed pump field is not chirped, we have T = 1/F. 

Coefficients A and B characterizing a coincidence- 
count pattern in a Hong-Ou-Mandel interferometer as 
given in Eq. (|40[) are determined for the general form of 
two-photon spectral amplitude $ as follows: 



A = 



B 



V(/2s+/2*/^-(/2^J2'''P 



(/is + fu? 



2 (/2s + /|i - f2si) ' 



2(/2s + /2*. + /2,s.) 

(A21) 
(A22) 



and coefficients Vfi- and £fr are given by formulas in 
Eqs. and 



APPENDIX B: SCHMIDT DECOMPOSITION OF 
TWO-PHOTON SPECTRAL AMPLITUDE $ 

In order to determine Schmidt decomposition of two- 
photon spectral amplitude $ as given in Eq. ([46]) we need 
reduced statistical operators belonging to the signal and 
idler fields. Statistical operator ps of the signal field can 
be written in the form: 



duj'^ / dujs'^s{i^s,^^s)al{ujs)\va,c){va.c\as{uj. 



(Bl) 

Weighting function determined along the formula 

nOO 

*,(w;,w^)=/ duj,<S>{uj',,uj,)<S>*{uJs,uj.,) (B2) 



takes the following form using Eq. ([9]) for the two-photon 
spectral amplitude $: 



X exp {eiLOs + el^^'s) 



(B3) 



and 



62 



e2c 



ei 



8/2" 



l/2siP 



/i 



flif2si 



(B4) 



Normalization constant introduced in Eq. (jB3[) . 

(el? 



1^ |2^./2(e^-e2. 



■ exp 



2(e^ - e2c 



(B5) 



guarantees normalization of the signal-field statistical op- 
erator Ps such that J_ dujs'^(ujg,LUg) — 1. Statistical 
operator pi of the idler field can be expressed similarly 
as that for the signal field. 

Coefficients An and functions (f>s,n and occurring 
in the Schmidt decomposition in Eq. (j46p are given as 
solutions of the following integral equations: 



duJa^a{uj'ai^a)(f>a,n{i^a) = A^(/)a^„(a;^), a = S,Z. 

_ J 

(B6) 

Using linear substitution the kernel in Eq. (|B3p can 
be transformed into the form: 

4'(x, y) = exp[-(l + P){x^ + y^) + 2xy] (B7) 



and 



P=^-l. 

e2c 



(B8) 
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It can be shown [46j that the following functions (/>„ 
obey the integral equation in Eq. (|B6|) for kernel de- 
fined in Eq. (|H7| : 



X H„ 



2"n!\/7n? 



exp 



1^ 
2d 



'l-?92 



1? 



n = 0,l,...,c»; (B9) 



symbols Hn denote Hermite polynomials. The corre- 
sponding eigenvalues form a geometric progression: 



0, 1, . . . , oo 



and parameter d is given as follows: 



d = l + P - VP2^^2P. 



(BIO) 



(Bll) 



A maximally entangled state according to entropy Se 
of entanglement determined in Eq. (|49|) is reached if all 
eigenvalues A„ are equal, i.e. when P — > 0. To under- 
stand this condition, we express coefficient P defined in 
Eq. (jBSP in terms of coefficients /: 




Re{/2s/|.'J) 



l.(B12) 



It can be shown that coefficient P goes to zero if coeffi- 
cient Vj defined in Eq. (jAll[) goes to zero too. Provided 
that coefficients Qs^ Qi, and Gsi can be omitted, the con- 
dition I?^ = is fulfilled only in the limit Tp — > and 
Zp 0. For this reason, we investigate the behavior of 
coefficient Vf with respect to pump-pulse duration Tp, 
pump-beam transverse width Zp, and widths as and 
of frequency filters in the area around Tp = or Zp = 0. 

For fixed values of Zp, as, and ai and around Tp = 
the following equality holds: 



dVf 



* SI p 



(B13) 



Also fixing values of Tp, ag, and ai and being around 
Zp = we arrive at: 



dVf 



d{Zl) 



y2r2 

^ SI p 



af 



y2 



VpsVp^gs^. (B14) 



If (dno(wa)/rfWaL^=^o)w° < Ua for a = s,i then g^si = 
2v^2s52i and the derivatives in Eqs. ()B13P and (|B14|) 
are positive (op = is assumed). This means that if the 
value of pump-pulse duration Tp is small [lower than that 
given by the condition in Eq. (|53p ] the shorter the pump- 
pulse duration Tp the more the down-converted fields are 
entangled. Similarly it holds for sufficiently small values 
of pump-beam transverse width Zp that the narrower 
the pump-beam transverse width Zp the more the down- 
converted fields are entangled. 

We get the following expressions for derivatives of co- 
efficient T) f with respect to widths CTs and ai of frequency 
filters: 



dVf 



9(^2) 



y2^2 



y2 2-2 



g^ 



(B15) 



For our waveguide, > and Gi > and this 
means that the derivatives of coefficient Vf written in 
Eqs. (jBlSp are negative (op = is assumed). Thus the 
wider the frequency filters, the smaller the value of co- 
efficient Vf and the more the down-converted fields are 
entangled. 



[1] C.K. Hong, Z.Y. Ou, and L. Mandel, Phys. Rev. Lett. 

59, 2044 (1987). 
[2] L. Mandel, E. Wolf, Optical Coherence and Quantum Op- 
tics (Cambridge Univ. Press. Cambridge, 1995). 
[3] D. Bouwmeester, J.-W. Pan, K. Mattle, M. Eibl, H. We- 

infurter, and A. Zeilinger, Nature 390, 575 (1997). 
[4] F. De Martini, V. Mussi, and F. Bovine, Opt. Commun. 

179, 581 (2000). 
[5] J. Pei^ina, Z. Hradil, and B. Jurco, Quantum Optics and 

Fundamentals of Physics (Kluwer, Dordrecht, 1994). 
[6] F.A. Bovino, G. Castagnoli, A. Ekert, P. Horodecki, CM. 

Alves, and A.V. Sergienko, Phys. Rev. Lett. 95, 240407 

(2005). 

[7] D. Bouwmeester, J.-W. Pan, M. Daniell, H. Weinfurter, 



and A. Zeilinger, Phys. Rev. Lett. 82, 1345 (1999). 
[8] D. Brufi and N. Liitkenhaus, in Applicable Algebra in 

Engineenng, Communication and Computing Vol. 10 

(Springer, Berlin, 2000); p. 383. 
[9] A. Migdall, Physics Today 52 (1), 41 (1999). 
[10] P.G. Kwiat, E. Waks, A.G. White, I. Appelbaum, and 

P.H. Eberhard, Phys. Rev. A 60, R773 (1999). 
[11] Y. Nambu, K. Usami, Y. Tsuda, K. Matsumoto, and K. 

Nakamura, Phys. Rev. A 66, 033816 (2002). 
[12] C.E. Kuklewicz, M. Fiorentino, G. Messin, F.N.C. Wong, 

and J. Shapiro, Opt. Express 13, 127 (2005). 
[13] S. Carrasco, J. P. Torres, L. Torner, A.V. Sergienko, 

B.E.A. Saleh, and M.C. Teich, Phys. Rev. A 70, 043817 

(2004). 



14 



[14] H.G. de Chatellus, A.V. Sergicnko, B.E.A. Salch, M.C. 

Tcich, and G. Di Giuseppe, Opt. Express 14, 10060 [31 
(2006). 

S.E. Harris, Phys. Rev. Lett. 98, 063602 (2007). [32 
X. Li, P.L. Voss, J.E. Sharping, and P. Kumar, Phys. 
Rev. Lett. 94, 053601 (2005). [33 
J. Fulconis, O. Ahbart, W.J. Wadsworth, P. S. Russell, 
and J. G. Rarity, Opt. Express 13, 7572 (2005). [34; 
M. Centini, J. Pei^ina Jr., L. Sciscione, C. Sibilia, M. 
Scalora, M.J. Bloemer, and M. Bertolotti, Phys. Rev. A [35 
72, 033806 (2005). 

J. Pefina Jr., M. Centini, C. Sibiha, M. Bertolotti, and 
M. Scalora, Phys. Rev. A 73, 033823 (2006). [36 
V. Giovannetti, L. Maccone, J.H. Shapiro, and F.N.C. 
Wong, Phys. Rev. Lett. 88, 183602 (2002). [37 
V. Giovannetti, L. Maccone, J.H. Shapiro, and F.N.C. 
Wong, Phys. Rev. A 66, 043813 (2002). [38 
O. Kuzucu, M. Fiorentino, M.A. Albota, F.N.C. Wong, 
and F.X. Kaertner, Phys. Rev. Lett. 94, 083601 (2005). [39 
V. Giovannetti, S. Lloyd, L. Maccone, and F.N.C. Wong, 
Phys. Rev. Lett. 87, 117902 (2001). 

A.B. U'Ren, K. Banaszek, and LA. Walmsley, Quantum [40 
Inf. Comp. 3, 480 (2003). 

A.B. U'Ren, C. Silberhorn, K. Banaszek, LA. Walmsley, [41 
R.K. Erdmann, W.P. Grice, and M.G. Raymer, Laser 
Phys. 15, 146 (2005). [42 
W.P. Gricc, A.B. U'Ren, and LA. Walmsley, Phys. Rev. 
A 64, 063815 (2001). [43 
S. Carrasco, A.V. Sergienko, B.E.A. Saleh, M.C. Teich, 
J.P. Torres, and L. Torner, Phys. Rev. A 73, 063802 [44 
(2006). 

J.P. Torres, F. Macia, S. Carrasco, and L. Torner, Opt. [45 
Lett. 30, 314 (2005). 

J.P. Torres, M. W. Mitchell, and M. Hendrych, Phys. Rev. [46 
A 71, 022320 (2005). 

A.B. U'Ren, R. Erdmann, and LA. Walmsley, J. Mod. 



Opt. 52, 2197 (2005). 

A.B. U'Ren, R.K. Erdmann, M. dc la Cruz-Guticrrcz, 
and LA. Walmsley, Phys. Rev. Lett. 97, 223602 (2006). 
J. Pefina Jr., M. Centini, C. Sibilia, M. Bertolotti, and 
M. Scalora, Phys. Rev. A 75, 013805 (2007). 
Y.J. Ding, S.J. Lee, and J.B. Khurgin, Phys. Rev. Lett. 
75, 429 (1995). 

A. De Rossi and V. Berger, Phys. Rev. Lett. 88, 043901 
(2002). 

M.C. Booth, M. Atature, G. Di Giuseppe, B.E.A. Salch, 
A.V. Sergienko, and M.C. Teich, Phys. Rev. A 66, 023815 
(2002). 

Z.D. Walton, M.C. Booth, A.V. Sergienko, B.E.A. Saleh, 

and M.C. Teich, Phys. Rev. A 67, 053810 (2003). 

M. Ravaro, Y. Seurin, S. Ducci, G. Leo, V. Berger, A. Dc 

Rossi, and G. Assanto, J. Appl. Phys. 98, 063103 (2005). 

L. Sciscione, M. Centini, C. Sibilia, M. Bertolotti, and 

M. Scalora, Phys. Rev. A 74, 013815 (2006). 

L. Lanco, S. Ducci, J.-P. Likforman, X. Marcadet, 

J.A.W. van Houwelingcn, H. Zbinden, G. Leo, and V. 

Berger, Phys. Rev. Lett. 97, 173901 (2006). 

Z.D. Walton, A.V. Sergienko, B.E.A. Saleh, and M.C. 

Teich, Phys. Rev. A 70, 052317 (2004). 

A.W. Snyder and J.D. Love, Optical Waveguide Theory, 

(Chapman & Hall, London, 1983). 

J. Pefina Jr., A.V. Sergienko, B.M. Jost, B.E.A. Saleh, 
M.C. Teich, Phys. Rev. A 59, 2359 (1999). 
C.K. Law, LA. Walmsley, and J.H. Eberly, Phys. Rev. 
Lett. 84, 5304 (2000). 

C.K. Law and J.H. Eberly, Phys. Rev. Lett. 92, 127903 
(2004). 

S. Parker, S. Bose, and M.B. Plenio, Phys. Rev. A 61, 
032305 (2000). 

V. Pefinova and A. Luks, private communication. 



